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We calculate the electromagnetic, axial and pseudo-scalar form factors of the Nucleon to 
A (1232) transition using two dynamical light degenerate quarks and a dynamical strange 
quark simulated with the domain wall fermion action. Results are obtained at lattice spac- 
ings a = 0.114 fm and a — 0.084 fm, with corresponding pion masses of 330 MeV and 
297 MeV, respectively. High statistics measurements are achieved by utilizing the coherent 
sink technique. The dominant electromagnetic dipole form factor, the axial form factors and 
the pseudo-scalar coupling are extracted to a good accuracy. This allows the investigation 
of the non-diagonal Goldberger-Treiman relation. Particular emphasis is given on the ex- 
traction of the sub-dominant electromagnetic quadrupole form factors and their ratio to the 
dominant dipole form factor, Rem and Rsm, measured in experiment. 

PACS numbers: 11. 15. Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70. Dj 

I. INTRODUCTION 

Form factors are fundamental quantities which probe the internal structure of the hadron. They 
are typically extracted from electromagnetic or weak scattering processes on hadronic targets, 
dominated by one-body exchange currents. The prime example are the form factors of the proton, 
which remain the most well-studied. Its electromagnetic (Sachs) form factors have been measured 
since the 50's pQ and static properties such as the magnetic moment and the charge radius are 
extracted. For recent reviews on the experimental and theoretical status we refer the reader to 
Refs. [U [2] and [3j H] respectively. Despite the long history of measurements of the electromagnetic 
nucleon form factors, polarization experiments recently revealed an unexpected behaviour in the 
momentum dependence of the electric to magnetic form factor of the proton which has triggered 
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theoretical investigations to explain the dynamics that give rise to such behavior [5]. 

The proton, being the building block of all matter that is presently observed to be stable, 
provides a nice laboratory for studying a relativistic bound state. One fundamental question is 
whether hadrons being composite systems are deformed and in particular whether the proton is 
spherical or has an intrinsic deformation. The elastic form factors do not suffice to answer this 
question on nucleon deformation, an important quantity that characterizes the distribution of 
quarks in the nucleon. The reason lies in the fact that the spectroscopic quadrupole moment of 
an J = 1/2 state vanishes identically in the laboratory frame if a one-photon exchange process 
is studied, although a quadrupole deformation may still exist in the body-fixed intrinsic frame. 
Therefore, regarding the nucleon, one has to study the transition to the lowest positive parity 
J = 3/2 state which is the A(1232). The jN A matrix element is parameterized in terms of a 
dominant magnetic dipole, Gmi, plus the sub-dominant electric quadrupole, Ge2, and Coulomb 
quadrupole, Gc2, transition form factors. Detection of non-zero Ge2 or Gq2 signals the existence 
of deformation in the N — A system [6HH] • Precise electroproduction experiments in the last decade 
demonstrated that this is indeed the case and provided measurements of the EM transition form 
factors for a wide range of values of the momentum transfer squared q 2 . The E2 and C2 amplitudes 
are measured to a few percent of the dominant, Ml, amplitude and are typically given as ratios to 
the Ml amplitude, denoted by Rem and Rsm respectively. 

State-of-the-art lattice QCD calculations can yield model independent results on hadron form 
factors, thereby providing direct comparison with experiment. Like in experiment, the electro- 
magnetic nucleon form factors have been studied by many collaborations recently using dynamical 
simulations [9Hl5j. Reproducing the experimental results on the electric and magnetic form fac- 
tors is a prerequisite for enabling lattice predictions of other form factors. This is also true for 
lattice calculations of the dominant magnetic dipole N to A transition form factor which is also 
well measured experimentally. In particular, in the case of the N to A, there are no disconnected 
contributions and therefore reproducing this form factor would provide a validation of lattice QCD 
techniques in calculating hadron form factors. The evaluation of the sub-dominant N to A elec- 
tric and Coulomb quadrupole form factors have also been studied for many years in dedicated 
experiments since, as we already pointed out, a non-zero value of these form factors signals a de- 
formation in the N-A system. However the experimental determination needs model input and 
therefore lattice QCD can provide an ab initio calculation of these fundamental quantities. 

In the axial sector, in the case of the nucleon, there exist two form factors, the axial, Ga, and 
induced pseudo-scalar, G p , form factors. They have been studied in neutrino scattering and muon 
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capture experiments, respectively but experimental data are less precise [TBI HZ] • There have also 
been several lattice evaluations of the nucleon axial charge gA [3 HBH2Q] and of the momentum 
dependence of the two form factors \12\ 121]. Partial conservation of axial symmetry (PCAC) leads 
to a relation between the nucleon axial charge and the pseudo-scalar ir — N coupling constant 
9nNN, the well-known Goldberger-Treiman relation. The strong decay of the A obscures greatly 
experimental studies of the N to A weak matrix element but some information on the dominant 
axial transition form factors C§(q 2 ) and Cg 1 (g ,2 ) is available from neutrino interactions on hydrogen 
and deuterium targets. C§ and C§ are the analogue of the nucleon axial form factors, Ga and 
G p , respectively. Indeed, like G p , the q 2 dependence of Cq is dominated by the pion pole and 
due to the axial Ward-Takahashi identity (AWI) a relation can be derived between C§ and the 
phenomenological strong coupling of the pion- nucleon- A vertex, g n NA- This relation is referred to 
as the non-diagonal Goldberger-Treiman relation. 

Such observations strongly motivate the study of the iV-to-A transition from first principles 
using lattice QCD. The first lattice study of the electromagnetic *yN A transition was carried out 
in the quenched approximation [22] at a fixed Euclidean momentum transfer squared Q 2 = —q 2 with 
inconclusive results as to whether the E2 or C2 amplitudes were non-zero due to large statistical 
errors. A study employing the formalism of Ref. [22J followed using quenched and two dynamical 
flavors of degenerate Wilson-type quarks at smaller quark masses but still only at the lowest q 2 - 
value allowed on the lattices at hand. Although there was an almost ten-fold increase in statistics 
the values obtained for the quadrupole form factors had large statistical noise and a zero value could 
not be excluded |23} 124] . In order to obtain sufficient accuracy we combined sequential inversions 
through the source instead of through the current for the evaluation of the three-point functions 
and optimized sources that led to a large sample of statistically independent measurements for a 
given q 2 -value. The calculation, carried out in the quenched approximation, confirmed a non-zero 
value with the correct sign for both of the quadrupole amplitudes [25\ |2"B] . A similar study was 
also carried out for the axial vector N to A matrix element [27]. Using this new methodology 
we extended the calculation of the N to A electro- weak form factors to unquenched lattice QCD. 
For the latter study we used Nj = 2 Wilson fermions as well as an Nf = 2 + 1 calculation with a 
mixed action with domain wall valence quarks on a staggered sea reaching a pion mass of about 350 
MeV [211 I28H30] . This calculation showed that the unquenched results on the Coulomb quadrupole 
form factor at low q 2 decreased towards the experimental results. However, the discrepancy in the 
momentum dependence of the dominant dipole form factor remained with lattice results having 
smaller values at low g 2 -values and a weaker dependence on q 2 . Using the same set of sequential 
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propogators as in the electromagnetic case the axial and pseudo-scalar N to A form factors were 
studied [SUED]. The strong coupling constant g n NA and non-diagonal Goldberger-Treiman relation 
were examined in detail and it was demonstrated that the behaviour is very similar manner to the 
corresponding relations in the nucleon system. 

In this work we study the iV-to-A transition using Np = 2 + 1 dynamical domain wall fermions 
simulated by the RBC-UKQCD collaborations [31]. This eliminates ambiguities about the cor- 
rectness of the continuum limit due to the rooting of the staggered sea quarks and the matching 
required in a mixed action. Preliminary results have been presented in Refs. [32]. We use two 
ensembles corresponding to lattice spacing a = 0.114 fm and a = 0.084 fm and physical vol- 
ume of (2.7 fm) 3 . Both lattice spacings are smaller than the lattice spacing used in our previous 
mixed-action calculation. This allows, for the first time, the investigation of cut-off effects on these 
hadronic observables. For each lattice spacing, we chose to perform the calculation on the lightest 
pion mass set available, namely at 330 MeV pions for the coarse lattice and 297 MeV for the fine 
one, in order to be as close as possible to the physical regime. The goal is, first, to check whether 
lattice results on the well measured experimentally dominant dipole form approach experiment. 
Secondly, we would like to see the onset of the large pion cloud contributions to the quadrupole 
form factors as predicted by chiral effective theory [33]. Thirdly we will extract the axial N to A 
coupling that enters in chiral expansions of the nucleon axial charge as well as the strong coupling 
constant g^NA- Determining these quantities together with the corresponding quantities gA and 
g n NN for the nucleon as well as for the A on the same gauge configurations will enable simultaneous 
chiral extrapolations to the physical point and yield more reliable results on these fundamental 
quantities. 

The paper is organized as follows: In Section [TT] we describe the general lattice setup and outline 
the techniques utilized to extract all the transition form factors from three-point functions measured 



on the lattice. In Section III we present in detail the decomposition of the electromagnetic N to 
A matrix element on the hadronic level in terms of the Sachs form factors and discuss the results 



for the electromagnetic transition form factors. In Section IV we give the corresponding matrix 
element for the electro-weak transition and discuss the results on the axial and pseudo-scalar form 
factors. Finally, the last section contains our conclusions and an outlook regarding further studies 
in the subject. 
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II. LATTICE SETUP AND TECHNIQUES 

We use the Nf = 2 + 1 dynamical domain wall fermion (DWF) ensembles generated by the 
RBC and UKQCD collaborations |3R 1541 [55] with the strange quark mass fixed at the physical 
point. Specifically, we consider gauge configurations on lattices of volume 24 3 x 64 corresponding 
to a pion mass of about 330 MeV and inverse lattice spacing a 1 = 1.73(3) GeV and 32 3 x 64 
corresponding to a pion mass of about 297 MeV and a' 1 = 2.34(3) GeV. We refer to the former 
lattice corresponding to a -1 = 1.73(3) GeV, as the coarse DWF lattice, and the one corresponding 
to a' 1 = 2.34(3) GeV, as the fine DWF lattice. 

Domain wall fermions preserve chiral symmetry in the infinite limit of the fifth dimension, 
L5. In actual computations L5 is finite leading to an additive contribution to the quark mass 
as defined through the Axial Ward-Takahashi Identity (AWI). For the coarse ensemble a residual 
quark mass of am res = 0.00315(2) has been measured by UKQCD-RBC [31] with the extent of the 
fifth dimension set to L5 = 16. The same L5 extent for the fine ensemble leads to a much smaller 
violation, measured to am res = 0.000665(3), or just 17% of the bare quark mass [TT] . 

Details about the lattice parameters used in this study are provided in Table |IJ where for 
comparison the relevant values of the parameters used in our previous study using the mixed 
action [12j [211 12H] are also given. 

Volume jV c d ° n 7 B (AWas.) N^t (^meas.) a" 1 [GeV] Z v Z A m u , d /m s [GeV] m N [GeV] m A [GeV] 

coarse Np = 2 + 1 DWF [31] 

24 3 x 64 200 (800) 398 (1592) 1.73(3) 0.7161(1) 0.7161(1) 0.005/0.04 0.329(1) 1.130(6) 1.457(11) 

fine Np = 2 + 1 DWF p] 

32 3 x 64 176 (704) 309 (1236) 2.34(3) 0.7468(39) 0.74521(2) 0.004/0.03 0.297(5) 1.127(9) 1.455(17) 

Hybrid action [12] 
DWF valence: am Ujd = 0.0138, am s = 0.081 

28 3 x 64 300 (300) 300 (300) 1.58(3) 0.01/0.05 0.353(2) 1.191(19) 1.533(27) 

TABLE I: Parameters for the calculation of the electromagnetic and axial transition form factors. The 
mixed action results from Refs. [3TJ[5H] are also included for completeness. In the second (third) column we 
show the number of gauge configurations used for the dominant (suppressed) form factors. For the DWF 
lattices the numbers in the parentheses next to the number of configurations are multiplied by four, since 
the coherent sink method was employed, showing the actual number of measurements taken into account 
in the overconstrained analysis. In the fifth and sixth columns we list the values of the vector and axial 
current renormalization constants, respectively that have been used as input parameters in our calculation, 
since we have used local currents and not the lattice conserved ones. 
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In order to create the proton and A + states we use the standard interpolating operators 

x P( x ) = e *bc [/«( x) C lb d b {x)] U C {x), (1) 



u Ta (x) C-f a d b (x)^ u c {x) + \u Ta {x) C 7(7 u fe (x)] (f(x)} 



(2) 



respectively. The J = 3/2 A state is described by the Rarita-Schwinger vector-spinor where 
a = 1, 2, 3, 4 is the Lorentz vector field index. C = 7472 is the charge-conjugation matrix. 

Form factors of the N — A transition are extracted on the lattice from the three-point function 

(G^ JV (t 2 ,t 1 ;p / ,p;r T ))= Yl e- l P'- X2 e + ^ Xi r^^|r[ X ^(x 2 ,t 2 )J At (x 1 , il )x^(O,0)]|fi) (3) 

X2, Xi 

In this notation, an initial nucleon state with momentum p is created at time zero and propagated 
to a later time t\ at which it couples to the current J causing a transition to the A state of 
momentum p' which is annihilated at a later time t 2 . q = p' — p is the momentum transfer. The 
projection matrices T T are given by 







ft 
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[ 1 


{0 





(4) 



The one-body currents considered in this work include the local vector current 

2 1 - 

V»( x ) = ^u(x)^u(x) - -d(x)j^d(x) , 

the axial-vector current and pseudo-scalar density 

A^x) = ip(x)^j 5 —ip{x) , P a {x) = ^(x)j 5 —ip(x) 



(5) 



(6) 



with r a the three Pauli-matrices acting in flavor space and tp the isospin doublet quark field. Note 
that due to the AJ = 1 nature of the transition, only the isovector part of contributes and, 
due to isospin symmetry, only the flavor diagonal operator r 3 needs to be evaluated. Inclusion of 
baryon states in the three-point function Q and the use of standard Euclidean spin-sums for the 
Rarita-Schwinger field 

-*7 • p + m/s. 



2mA 



and the Dirac spinor 



u(p, s)u(p, , 



2p a p T ,p a ^ T -p T ^ a 1 

S °T+ -=—2- - 1 O 07^7^ 

3mi 3mA 3 



-»7 • p + mjv 



2mjv 



(7) 



(8) 



7 



lead to the isolation of the desired matrix element, assuming that the initial and final ground states 
dominate the propagation before and after the operator insertion, respectively. In order to cancel, 
in the large Euclidean time limit, the dependence on the Euclidean time evolution and on the 
unknown overlaps of the nucleon and A states with the initial states, we form the following ratio: 



,,,A./„A 



(Gt A (t 2 ,p';r 4 )> 

(G£A( t2i p' ; r 4 )) (G^ih-t^u)) (G£A( tli p/ ; r 4)) i V2 

<G*W(i 2 ,p;r 4 )) (Gt A (t 2 -ti,p';r 4 )) (G^(t 1 ,p;U))_ 
which requires also measurements of the nucleon (G NN ) and A (G^r) two-point functions 



(9) 



(G NN (t,p;T)) =J2e- ipX (0| T X *(x, t)^(0, 0) |0 ), (10) 

X 

<G A T A (t, P; r)) = ^ e -P-x r ^ <Q|T X? (x,i)^(0,0) |n). (11) 

X 

Implicit summations on indices % = 1,2,3 are assumed in the above ratio ([9]), which is designed 
such that the time evolution (and consequently the noise) appearing in its two-point function part 
is minimized. In the large Euclidean time limit (i 2 — t\ 3> 1, ti S> 1) where we have ground state 
dominance this ratio ^ thus yields a a time- independent function H"^ (p ', p ; T T ; [a) that is related 
to the matrix element (A(p')| J\n(p)). Therefore we look for the plateau region of Eq. Q in order 
to extract the matrix element that we are interested in. For a given operator insertion J and 
projection matrix IV, the function II^p'jp ;T T ; fi) is a linear combination of the corresponding 
form factors. These relations for the appropriate choice of T T and A vector index a will be given 
in the following sections. 

The computationally intensive part of the calculation lies in the calculation of the three-point 
function given in Eq. ^ . In order to achieve the extraction of the momentum dependence of the 
matrix element for the V^(x), A^(x) and P 3 (x) insertions, one needs an evaluation for a large 
number of values of the momentum transfer q. This is feasible by evaluating the matrix element 
using sequential inversions through the sink. In this method, the quantum numbers of the source 
and sink interpolating fields are fixed, effectively by fixing the a and r indices. The time slices of 
the source and sink are, in addition, fixed. The quark propagator with the operator insertion is 
obtained by the joining of a forward propagator and the sequential propagator which is obtained 
by using as a source the baryon state at the sink folded in with the two forward propagators 
from the source. With the forward and sequential propagators available, the operator insertion at 
selected intermediate times t\ and momenta transfers q is readily available. In this method the 
final state, in this case the A-state, is always at rest. Since the a — r space of indices still spans 
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a set of 16 independent inversions that would be required, an optimization in this space has been 
exploited. Three linear combinations are constructed from which the EM, axial and pseudo-scalar 
form factors are extracted such that the maximal set of statistically independent measurements of 
momentum transfer vectors q per q 2 value is achieved. In addition, they are chosen to decouple 
the dominant dipole (Ml) part of the EM transition from the sub-dominant quadrupoles E2 and 
C2 measurements. The three linear combinations which we construct and measure in this work 
are given below. 

3 

5 1 J (q;J)=^n^(0,-q;r 4 ;J) (12) 

<T=1 

3 

si(^,J) = E n^(o,-q ; r fc ;j) (13) 

<7^fc = l 

Si(q; J) = U J 3 (0, -q ; r 3 ; J) - \ [llf (0, -q ; T i; J) + n 2 >, -q ; T 2 ; J) J , (14) 

where J denotes the operators V^, and P 3 . Occasionally we refer to Si, S^, £3 a s optimal A 
sinks, although they actually correspond to an optimal linear combination of the full N — A three- 
point function with arbitrary insertion J. We stress that, given the forward propagators, three 
inversions in total are required in order to compute the momentum dependence of the full N — A 
transition and extract the electromagnetic, axial and pseudo-scalar form factors. 

Since the source-sink separation is fixed in this method it is crucial to suppress the excited baryon 
states as much as possible. This is achieved by employing gauge invariant Gaussian smearing on 
the local quark fields with APE-smeared gauge fields and parameters that have been carefully 
optimized for the nucleon state. For the coarse lattice, we show in Fig. [I] a comparison of results 
obtained with a sink-source separation of 0.91 fm and 1.14 fm. As can be seen, extending the 
source-sink separation to 1.14 fm, the plateau values for the dominant magnetic dipole form factor 
Gmi, which are the most accurate, are consistent with a time-separation of 0.91 fm. Since the 
larger time separation introduces a doubling in the statistical noise, for the accuracy needed in 
this study, we opt to use the smaller sink-source separation in time. For the fine lattice, we take a 
sink-source separation of AT = 12a corresponding to 1.01 fm, consistent with our findings using 
the coarse lattice. 

In order to improve accuracy, a goal that is particularly crucial for the extraction of the sub- 
dominant electromagnetic form factors, we employ a new method first implemented in the study of 
the nucleon form factors [11] and referred to as the coherent sink technique. The method consists 
of creating four sets of forward propagators for each configuration at source positions separated in 
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FIG. 1: The ratio R% from the source Si of Eq. ^ versus t/a for a source-sink separation 0.91 fm shifted 
by a time slice (triangles) and 1.14 fm (squares) for the four smallest non-zero q 1 values. The fit range is 
also shown along with the fitted lines and the corresponding error bands. The behavior is the same for both, 
but the error reduction is better in the former, which is what we therefore utilize in the calculations. 



time by one-quarter of the total temporal size. Namely, for the coarse DWF lattice, Nl = 24, we 
have forward propagators generated with sources positioned at: 

L L 3L 

(0,0), (-,16a), (-,32a), ( T ,48a) 

and for the fine DWF lattice, Nl = 32, placed at: 

(0,10a), (^,26a), (0,42a), (~,58a)}. 



or 



(-,10a), (~ 26a), (-,42a), (^,58a)|. 

From each source (xi,Ti), a zero-momentum projected A source is constructed at To slices 
away, i.e. at (xj,T + To). For the coarse DWF lattice To/a = 8, while for the fine DWF lattice 
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To I a = 12. Then a single coherent backward propagator is calculated in the simultaneous presence 
of all four sources. The cross terms that arise vanish due to gauge invariance when averaged over the 
ensemble. The forward propagators have already been computed by the LHPC collaboration [TT] 
and therefore we effectively obtain four measurements at the cost of one sequential inversion. This 
assumes large enough time-separation between the four sources to suppress contamination among 
them. A question that arises is whether or not there exist statistically important correlations 
among these four measurements. In Fig. [2] we show the dependence of the jackknife error on the 
magnetic dipole Gmi for different coherent sink bin sizes. As can be seen, the jackknife errors 
using one sequential inversion for each are the same as combining all four in single inversion. This 
is a direct verification that cross-correlations between the different sinks are absent or negligible. 
2> 1 — 



1.5 - 



0.5 - 



T 



T 



T 



o Coherent sink data treated independently, 800 data points 
* 2 coherent sink points averaged, 400 data points 
» All coherent sink points averaged, 200 data points 



H 



i 







I 



4 



Q 2 (GeV 2 ) 



FIG. 2: Dependence of the jackknife error for Gmi{Q 2 ) on the coherent sink bin sizes. This test shows that 
there is no problem with cross-correlations in the coherent sink method applied in this study. 



Finally the full set of lattice data obtained at a given Q 2 value is analyzed simultaneously 
by a global x 2 minimization using the singular value decomposition of an overconstrained linear 
system |21| 136] . Generically, this consists of setting up the following linear over-complete system 
of equations 

P(q; fJ ,)=D(< l ; f j,)-F(Q 2 ), (15) 
where -P(q; n) represent the lattice measurements of the appropriately defined ratios of Eq. [91 each 
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one with its associated statistical weight Wk- The column vector F(Q 2 ) contains the number M 
of form factors to be extracted. If we let N represent the number of momentum vectors q and 
current directions \i that contribute to a specific value of Q 2 , then D(q; n) is a matrix structure of 
the form N x M which depends on kinematical form factors obtained from the trace algebra on 
the employed matrix element. The form factors, at the specific Q 2 value, are then extracted from 
the minimization of the total x 2 : 




by applying the singular value decomposition on the N x M, D(q; n) matrix. All the errors on the 
lattice measurements as well as the errors on the form factors are determined from the jackknife 
procedure. 



III. ELECTROMAGNETIC N TO A TRANSITION FORM FACTORS 
A. The electromagnetic matrix element 

The electromagnetic transition matrix element 

{A( P ',s')\j,\N( P ,s)) =i^| ^-^^^y\ a {p',s')O a ,u(p,s) (17) 
is decomposed in terms of three multipole form factors: 

O afl = G M1 (q 2 )K^ + G E2 (q 2 )K* 2 + Gciiq 2 )^ 2 

where the kinematical factors in Euclidean space are given by 

k mi 3 m A + m N lP 

K ^ ~ (m A + mN )2 + Q2 2m N 1£ ^ p p ' 

K% = ~K% + 6 n-' (Q 2 ) "^^2 i 75 e aXa , p^e^ p,p> 5 , 

K% = -6""V) mA 2 *™" i75 fb (Q 2 (P + P% ~ q ■ (p + p'K) ■ (18) 

Thep(s) andp'(s') denote initial and final momenta (spins), q 2 = (p' —p) 2 , and u a (p', s') is a Rarita- 
Schwinger vector-spinor. We also define Q(Q 2 ) = [(ttoa + rn^) 2 + Q 2 ] [(ttt-a — ttin) 2 + Q 2 ] , with 
(Q = q, Q4 = iq°), so the lattice momentum transfer gives Q 2 = —q 2 . 

In this work we present results for the dominant magnetic dipole form factor Gmi{q 2 ) as well 
as the sub-dominant electric Ge2(q 2 ) and Coulomb quadrupole Gc2(q 2 ) form factors. Note that 
these are all scalar functions depending on the momentum transfer q 2 = —Q 2 , whereas on the 
lattice only the space-like q 2 are accessible, thus Q 2 > 0. 



12 



B. The magnetic dipole form factor 

The magnetic dipole form factor is directly evaluated from the optimized linear combination 
SY with the vector current V^x) insertion. In the large Euclidean time separation limit with the 
A produced at zero momentum we obtain, 

SXteVp) = »a|(p2-P3)^ + (P3 -Pi)h,i* + (pi -P2)^G M i(Q 2 ) ■ (19) 

The vector index \i takes spatial values, /i = 1, 2, 3 and A is a kinematical constant, 

m A + m N J—Ejj—^ 
AitinEn V En + tun 

The local vector current of Eq. ([5]) is not conserved by the lattice action and the renormalization 
constant Zy, given in Table [IJ is used to renormalize the current. Zy is determined from charge 
conservation that dictates that the electric nucleon form factor is one at Q 2 = 0, namely Zy = 
1/Ge(0) = 1/-Fi(0) where F\ is the Dirac form factor. 

In Fig.[3]we show the DWF results for the magnetic dipole form factor Gmi at m n = 330 MeV on 
the coarse lattice and at m n = 297 MeV on the fine lattice as a function of the momentum transfer 
Q 2 . These are compared with our previous results obtained with a hybrid action approach that 
uses Asqtad improved staggered fermions generated by the MILC collaboration |37j and domain 
wall valence quarks [28] . The experimentally available data (for more details see Ref. [28] ) are also 
shown in Fig. [3] showing a discrepancy between lattice results and experiment. Although there is a 
small decrease in the value of Gmi at high Q 2 bringing lattice data closer to experiment, the slope 
at small Q 2 is still smaller than in experiment. The second observation is that although the hybrid 
calculation used a lattice spacing about 50% larger than the fine DWF lattice, these data show no 
significant finite a-effects. Fits to a dipole form, <?o/(l + Q 2 / m o) 2 ' as wen as t° an exponential form 
go exp(— Q 2 /rho) are shown for the fine DWF lattice. As can be seen they both provide a good 
description of the lattice results. A list of the fit parameters for all sets is provided in Table |Tlj 

The discrepancy between experiment and lattice results is clearly reflected in the value of the 
dipole mass of mo = 0.78 GeV obtained by performing a dipole fit to the experimental data as 
compared to the mo values obtained from the lattice results listed in Table [TTJ The steeper rise of 
the experimental results on Gmi as a function of Q 2 near the origin is indicative of the onset of 
strong chiral quark effects, or equivalently, the lack of strong pion cloud from the still heavy pion 
mass lattice ensembles that are utilized. Similar behavior has also been observed in the nucleon 
electromagnetic form factors studies |10[ [TT| 113] , The N to A transition is particularly clean 
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FIG. 3: The magnetic dipole Gmi(Q 2 ) using DWF fcrmions (both coarse and fine lattices) and using the 
hybrid action. The circles show the experimental results. The solid blue (dashed black) line is a fit to dipole 
(exponential) form for the fine DWF lattice. 



since there is no ambiguity regarding disconnected contributions and thus the natter dependence 
observed in the N-to-A electromagnetic form factor must be of different origin. Large pion cloud 
effects would have to set in as we lower the pion mass in order to explain the experimental curve. 
Such effects have been shown to arise in chiral expansions |33j and it is thus interesting to repeat 
the calculation for < 250 MeV where they are expected to become more pronounced. 
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m n [GeV] 


g„ m [GeV] 


9o 


too [GeV] 


coarse Np 


= 2 + 1 DWF 






0.329(1) 


1.937(65) 1.171(44) 


1.737(53) 


1.025(32) 


fine Np = 


2 + 1 DWF 






0.297(5) 


2.115(161) 1.078(79) 


1.907(127) 


0.939(55) 


HYBRID 








0.353(2) 


3.263(64) 1.305(27) 3.05935(61) 1.097(21) 


Experiment- (various sources) 








3.266(40) 0.745(4) 


2.202(67) 


0.776(10) 



TABLE II: The fit parameters for the magnetic dipole form factor obtained for both a dipole and an 
exponential ht form with fit parameters (go, too) and (go, mo), respectively. 

C. The electric quadrupole form factor— Gei 



The sub-dominant electromagnetic quadrupole form factors Ge2 and Gc2 are extracted from 
the optimized sources and The relevant expressions for a static A final state are [28]: 



^(q;^) = - 3^< (0? 2 +Pa)8i,p + (P3+Pi)h,^ + {pi+P2)h, l ?)QE2{Q 2 ) 



T) 

- 2-|(piP2 + PlP3 + P2P3} 



Ge2(Q 2 ) + ^f^ 2 (Q 2 ) 



(21) 



for the spatial current directions fj, = 1, 2, 3. For the temporal current direction fi = 4, we have 



^2 / (q;/^ = 4) 



-i 6 B 



(PlP2 + PlP3 + P2P3)QC2{Q 2 ) , 



(22) 
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where B is given by B = g^— -A, and A is the constant provided in Eq. 

Notice that the above combination, if used alone, will not allow for the extraction of Gc2 at the 
lowest photon momentum q = (1,0, 0)^. Since chiral effects are stronger at low Q 2 values and 
experiments are targeted in that regime, we utilize the optimal linear three-function combination 
SY m order to obtain Gc2 also at the lowest Q 2 point allowed on the lattice. The corresponding 
expressions are 



Safa P) 



3 A 



Pn 



3 5. 



pI 



2mA 



1-3 4i <j<-2(cr 



(23) 



for \i = 1, 2, 3 and for the temporal component 

^(q;M = 4) = ^ 



1 - 3 



pi 



Gc2(Q 2 ), 



(24) 
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which is directly proportional to Gc2{Q 2 )- Data obtained from both S% and S% are simultaneously 
fitted in the overconstrained analysis in order to extract the momentum dependence of Ge2 and 
Gc2 as accurately as possible. 
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FIG. 4: In (a) the result of the electric quadrupole form factor Gei{Q ) extracted from the fine DWF lattice 



measurements is shown. The results obtained from the hybrid action 128], as well as the experimentally 
extracted results from Bates [381(40] , Jlab [41] and MAMI [42] [43] are also plotted for comparison. In |(b) 
the corresponding Rem evaluated in the rest frame of the A baryon (p' = 0) is depicted for the fine DWF 
lattice as well as for the hybrid action [55]. The experimentally available results from [58H13"] are also shown. 



In Fig. |4(a)| we plot the values of the electric quadrupole form factor Ge2 for a range of values 
of Q 2 < 1 GeV 2 , in the case of the fine DWF lattice. These results are compared to the results 
obtained from the mixed action |28j. We also mention here that in the case of the coarse DWF 
lattice the statistical noise on the Ge2 and Gc2 values is large, so a zero value can therefore not 
be excluded. The phenomenologically interesting ratio Rem is defined as 

G E 2{Q 2 ) 



Rem 



(25) 



Ga/i(Q 2 )' 

and has been used traditionally as a signal of deviation from spherical symmetry in the nucleon-A 
system. Early quark models as well as models of the proton wave function based on relativistic 
quarks including two-body exchange currents agree that a small Rem value in the — 1 ~ 2% regime 
should appear. The experimental values included in Fig. |4(a)| show practically no dependence on 
Q 2 . The same is true for the lattice data and in fact a good consistency with the experiment is 
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evident. The approach to the physical point can be predicted in chiral effective theory |33j where 
a non-monotonic dependence on the pion mass is expected with a minimum at 200 MeV. It is a 
significant challenge for the lattice to provide accurate results in the future in this regime in order 
to crosscheck the pion dynamics. 



D. The Coulomb quadrupole form factor— Gci 

As mentioned in the previous section, the Coulomb quadrupole form factor is computed with 



the help of Eqs. (21)-(24). In the case of Gc2, Fig- 5(a) shows the results from the fine DWF 
lattice for values of Q 2 < 1.5 GeV 2 . The values of Gc2(Q 2 ) are positive and consistent with 
previous results obtained using the mixed action [28J, and are also shown on the same figure. The 
experimentally measured ratio of Coulomb quadrupole to magnetic dipole form factor known also 
as CMR is defined by 

„ _ |q| Gc2(Q 2 ) . . 

in the frame where the A is produced at rest. Lattice results on the Rsm ratio are shown in Fig. |5(b)"] 
where tba in Eq. (|26|) is set to the physical mass. Known values of Rsm from various experiments 



are included in Fig. 5(b) and as with Rem show almost no dependence on the momentum transfer. 
This is also the feature shown by the two lattice ensembles, the fine DWF at 297 MeV and the 
hybrid scheme at 353 MeV which are in very good agreement with each other. Despite the large 
statistical errors which escort the lattice values, they disagree with the experiment. Chiral effective 
theory predicts a monotonic decrease of this ratio as the pion mass approaches the chiral limit, 
which is different from the dependence of Rem- The onset of large pion effects are expected below 
300 MeV pions. 

The overall conclusion is that QCD confirms non-zero quadrupole amplitudes pointing to the 
existence of the deformation in the N — A system, as coded in the EMR and CMR ratios. However, 
quantitative agreement with experiment has to await simulations at lighter pions masses, expected 
to become available in the next couple of years. Statistical accuracy at these light pions masses in 
particular for the sub-dominant form factors is an issue that has to be addressed. The use of the 
coherent source technique as employed here is a way to increase statistical accuracy. 
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FIG. 5: In plot (a) we show the Coulomb quadrupole form factor Gc2(Q 2 ) extracted from the fine DWF 
lattice measurements. Along with it we provide also the result from the hybrid action approach [28 . Plot |(b) 
depicts the corresponding Rsm evaluated in the rest frame of the A baryon. Non-zero values are confirmed, 
for the lowest Q 2 values accessible on the lattices. We also show results using the hybrid action taken from 
Ref. . Experimental results are also included using the same notation as those in Fig. |4j 



IV. AXIAL N TO A TRANSITION FORM FACTORS AND THE 
GOLDBERGER-TREIMAN RELATION 

A. The Electro-weak and Pseudo-scalar transition matrix element 



The nucleon to A matrix element of the axial vector current is parameterized in terms of four 
dimensionless form factors. In the Adler parameterization |44j it is written as follows 

1/2 



(A(p',s')\Al\N(p,s)) 

T • 9 " 

mi 



m&rriN 



3 V£a(p')£aKp) 



m 



N 



u P (p, s) (27) 



TUN no N 

with the axial current given in Eq. ([6]). 

The form factors C^(q 2 ) and C^(q 2 ) belong to the transverse part of the axial current and are 
both suppressed |27j relative to the longitudinal form factors C§(q 2 ) and C§(q 2 ), which are the 
dominant ones and are the ones considered in this work. 
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Likewise, the pseudo-scalar transition form factor GWa(<7 2 )) is denned via 



2m g {A(p',s')\P 3 \N{p,s)) = ( 



3 \E A (p')E N (p) 



1/2 



2 _ 2 U A+(P > s )7^T u P\Pi s ) 



rat 



2m^ 



(28) 



where the normalization of the RHS of (28) is chosen such that G n NA{q ) reproduces the phe- 



nomenological coupling of the it — N — A vertex in the strong interaction Lagrangian, 

9-kNA 



NA 



_ } -A^TT -T N + h.C. 



(29) 



and the pseudo-scalar density is defined in Eq. ([6j). In the SU(2) symmetric limit with m q denoting 
the up/down mass, the pseudo-scalar density is related to the divergence of the axial- vector current 
through the axial Ward-Takahashi identity (AWI) 



0M£ = 2m q P a 



(30) 



Taking matrix elements of the above identity between N and A states leads to the non-diagonal 
Goldberger-Treiman (GT) relation 



c 5 V) + \c£{f) 



1 G nN A(q 2 )Uml 



m% 



mt 



(31) 



o N 2m N 

On the other hand, flavor symmetry in the hadronic world is expressed through the partially- 
conserved axial vector current (PCAC) hypothesis 



2_a 



0M£ = f n mi7T 



(32) 



which relates the pseudo-scalar current to the pion field operator and the pion decay constant 



which is here is taken to be 92 MeV. From Eqs. (30) and (32) the pion field 7r a is related to the 
pseudo-scalar density via 



7T 



2m q P a 



Assuming pion pole dominance we can relate the form factor C§ to G w na through: 

_}_r<Ai 2\ 1 G n NA(q 2 )fn 
tun 2 m% — q z 



(33) 



(34) 



Then, substituting Eq. (34) in Eq. (31 ), we obtain the simplified Goldberger-Treiman (GT) relation 



GWa (q Z )U = 2m N C£(q 2 " 



(35) 
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in an analogous fashion to the well known GT relation which holds in the nucleon sector studied on 
the lattice in Ref. [21]. Pion pole dominance therefore fixes completely the ratio C^(g 2 )/C^(g 2 ) 
as a pure monopole term 



in 



N 



m 2 — q 2 



(36) 



The aim here is to calculate the dominant axial C§(q 2 ), Cq(Q 2 ), as well as the pseudo-scalar 
G^NAiQ 2 ) form factor and examine the validity of the GT relations within the dynamical DWF 
framework, using both the coarse and fine DWF lattices. 



B. The dominant axial C§, C$ transition form factors 



The extraction of the axial transition form factors requires data from two sets of the optimal 
A sinks, namely S± and S2, which are introduced in section II, for the local isovector axial-vector 
current insertion A^(x). The corresponding expressions for the large Euclidean time separation 
ratios are: 



Sf(q;j) =b 



{{E N -2m A + m N )+ (j^ 
^ \k=i 



P 



— (E N - m A )C£ + m N C£ - °^ 
rriN m,N 



for spatial components j = 1, 2, 3 of the axial current, and 

3 



En + rn-N 



(37) 



P 



k=l 



m N 



m N 



(38) 



for the temporal component. Since the four form factors are not completely decoupled by the above 
relations, we also employ the optimal A sink given in the plateau by 

3 \ 



3A 



£p* (Mp 2 - p") + <^(p 3 - p 1 ) + s^ip 1 - p 2 )) a 



\k=i 



(39) 



valid for spatial components j = 1,2,3. The kinematical factors A and B are given by 

B 



A 



B 



2 y/(E N + m N ) /E N 



(40) 



(En + mjy) 7 V 3 3mjv 

Data from and determine all four form factors C3 , C4 , C§ and Cq at each value of Q 2 in 
a simultaneous overconstrained analysis. Za is required to renormalize the axial vector operator. 
This has been computed by the UKQCD-RBC and LHP collaborations for both ensembles |11| 131} 



20 



. The values provided in Table [T] confirm that Zy = Za + 0(a 2 ) in the chiral limit, as expected 
for the manifestly chiral DWF action. 

The results for the axial dominant form factor C§ from the two DWF lattices considered in this 



work are presented in Fig. 6(a) and are in good agreement with the results obtained from the mixed 
action approach at = 353 MeV [21]. The Q 2 dependence is well described by two-parameter 
dipole (solid line) and exponential (dashed line) forms do / (1 + Q 2 / m A ) 2 , do exp(— Q /rh A ), respec- 
tively, which are almost indistinguishable in the plot. The fitted values for C§(0) = do (or do of 



the exponential form) and the corresponding axial mass niA (wia) are given in Table III In the 
same figure, we also show a dipole fit to the available experimental data [IS] which determine an 
axial mass within the range of values of tua ~ 0.85 — 1.1 GeV [45 1 I46j . obtained from the pure 
dipole parameterization. As in the case of Gmi{Q 2 ), we observe a flatter slope for the lattice data, 
reflected in the larger value of the axial mass tua extracted from the lattice results. 

m w [GeV] m,A [GeV] do fh>A [GeV] do m [GeV] co A' SirJVA 

coarse Np = 2 + 1 DWF (K) (a') 

0.329(1) 1.588(70) 0.970(30) 1.262(36) 0.940(21) 0.509(15) 5.132(204) 0.030(5) 9.525(168) 13.936(588) 
fine Np = 2 + 1 DWF 

0.297(5) 1.699(170) 0.944(58) 1.314(98) 0.927(46) 0.507(33) 5.756(516) 0.037(6) 8.444(491) 16.257(867) 
Hybrid action 

0.353(3) 1.795(40) 0.903(11) 1.386(18) 0.888(8) 0.496(10) 5.613(150) 0.019(11) 9.323(219) 11.446(617) 

TABLE III: The first column gives the pion mass in GeV. The second and third columns provide the dipole 
fit parameters m A and d extracted from fitting to do/(l + Q 2 /m A ) 2 , the fourth and fifth columns the 
corresponding parameters obtained from the use of an exponential ansatz do exp(— Q 2 /ih A ), the sixth and 
seventh columns the fit parameters m and Co extracted from fitting the ratio Cq/C^ to a monopole form 
cq/(1 + Q 2 /m 2 ) for the N-to-A process. The eighth and ninth columns show the calculated values of the fit 



parameters a' and A' defined in the linear fit of Eq. (51 1. The last two columns give the predicted values of 



the strong coupling constant g n NA = G^ata(O). The first value of the strong coupling constant is determined 



using the fit function of Eq. (501, while the second uses the linear fit based on Eq. (511, which is exactly 
equal to a'. 



In Fig. 6(b) we show the ratio Cq/C§, since ratios of this type are expected to be less prone 



to lattice artifacts. The dashed black line shows the pion pole dominance prediction of Eq. (36), 
where for mjv and i% we use the lattice extracted values that correspond to the fine DWF lattice. 
The predicted curve does not describe the data at low-Q 2 i.e., in the regime where the strong pion 
cloud effects are expected to be present. However, the fit to the monopole form cq/(1 + Q 2 /m 2 ) 
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FIG. 6: Plot (a) shows the Q 2 -dependence of the axial form factor C§ extracted from the coarse and fine 
DWF lattices. The corresponding mixed action results [5T] have also been included. The solid blue (dashed 
black) line is from the dipole (exponential) fit for to the fine DWF lattice results. Note that the error band 
corresponds to the dipole fit. The dotted brown line is the dipole fit to the experimental data. The ratio 



Cq/C§ versus Q 2 is plotted in (b) The dashed black line refers to the fine DWF lattice results and is the 
pion pole dominance prediction of Eq. (36). The solid blue line is a fit to a monopole form co/(l + Q 2 /m 2 ). 



describes satisfactorily the ratio yielding a heavier mass parameter m than the lattice value of 



the pion mass (see Table III). Such behavior has been observed also for the hybrid and quenched 
Wilson actions 1211. 



The lattice results for the C§ are plotted on Fig. [TJ The curve shown (solid line) in the figure 
corresponds to the form 

do c 



(l + Q 2 A4) 2 (l + Q 2 /m 2 )' 



(41) 



where cq and m are the parameters of the monopole term given in Eq. (36) that are expected to 



describe well the Cq/C^ ratio provided the pion pole dominance is applicable. The form described 



by the expression of Eq. (41), seems to provide the best fit to the fine DWF data. On the other 



hand, Cg is related to the C§ form factor through the expression 



Qs (Q 2 ) = ct{Q 2 



rn 



N 



ml + Q 2 ' 

The curve that corresponds to the dashed line is obtained from fitting the fine DWF data to this 
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form. In this case C§ is being described by the dipole form shown in Fig. 6(a) , while the nucleon 
and pion masses are the lattice evaluated ones. 
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FIG. 7: Lattice results for Cg are shown as a function of Q 2 . The solid blue line is the fit to the form of 



Eq. ( 41 1 , while the dashed black line corresponds to the form ( 



ml+Q 2 



Note that for the latter fit, the 



C5 factor is described by the dipole fit parameters. 



C. The Pseudo-scalar transition form factor and Goldberger-Treiman relation 



The pseudo-scalar form factor G^na(Q ), defined via the matrix element given in Eq. (28), 



is extracted directly from the optimized linear combination S\ with the pseudo-scalar current 
operator insertion of Eq. Q. In the large Euclidean time limit where only the nucleon and A 
states dominate the corresponding ratio yields 



Sf(q ; 75) 



2 E N + rriN 
3V En 



q\ + 92 + <?3 



(42) 



6m n 2m q (m1 + Q 2 ) 

Notice that the extraction of G n NA from the above equation requires knowledge of the quark 
mass m q and the pion decay constant, /„-, on the given ensembles. Calculation of /„• requires the 
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two-point functions of the axial-vector current A\ with local-smeared (LS) and smeared-smeared 
(SS) quark sources, 

c£ s (t) = E (°l T ( A K*> ^4(0, o)) |n > (43) 

X 

(and similarly for Cg S ), where A|(x, t) denotes the local operator and A^(x,t) the smeared oper- 
ator. The pion-to-vacuum matrix element 

(Q\A«(0)\ir b (p))=iUp,6 ab (44) 

is extracted from the two-point functions C£ s and Cg S and 



ff(t) = Z M P^ ° Ls{t) e m ^ 2 . (45) 



yields f n in the large Euclidean time limit. 

The renormalized quark mass m q is determined from the AWI, via two-point functions of the 
pseudoscalar density with either local (P 3 ) or smeared (P 3 ) quark fields, 



cf s (t) = E <n| t (p 3 (x, t)F 3 (o, o)j |n ) , (46) 

X 

(and similarly for Cg S ). The effective quark mass is defined by 



aw i ...rrhrZA C£ s (t) / Cf s (i) 



"*•« (t) --z;cEw\/c&w- (47) 

and its plateau value yields m g . Note that Zp will be needed only if ones wants m q alone. Since Zp 



enters also Eq. ( 42 ) it cancels -as does Za since it comes with f % - and therefore G w na is extracted 



directly from ratios of lattice three- and two-point functions without prior knowledge of either Za 



or Zp. We also note that the quark mass computed through (47) includes the effects of residual 
chiral symmetry breaking from the finite extent L5 of the fifth dimension. These effects are of the 
order of 60% for the coarse ensemble and 17% for the fine ensemble. Chiral symmetry breaking 



affects the PCAC relations and therefore the value of G-^na through Eq. (42) 
The ratio 

(Q 2 ) 

2m N C£(Q*) 



(48) 



is depicted in Fig. |8(a) It should be unity if the off-diagonal Goldberger-Treiman relation of 
Eq. (35 ) is satisfied, which in turn requires that PCAC holds exactly at the pion masses simulated in 
these ensembles. Deviations from this relation are seen in the low-Q 2 regime. For the fine ensemble 
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FIG. 8: In|(a)|we plot the ratio of Eq. ( 48 ) as a function of Q as a validity test of the GT relation. Similarly, 



in plot |(b)| ratio of Eq. ( 49 1 that relates to the validity of Eq. ( 34 ) 



considered in this study, the deviations from unity are less severe. At momentum transfers, of about 
Q 2 ^ 0-5 GeV 2 , the relation is at least approximately satisfied and it is consistent among all actions 
considered here. 



Pion pole dominance relates C§ to C§ through Eq. (36), which is very well satisfied by the 



lattice data for all the actions employed in this work (see also Fig. |6(b)D - This agreement is also 
evident from Fig. |8(b)[ where the ratio 

2(m2+Q2 )q A (g2) I ) 

is consistent with unity. 

In Fig. [9] we compare results on G^NAio 2 ) using the dynamical DWF lattices to the results 
obtained from the hybrid scheme taken from Ref. [3T]. There is an agreement for Q 2 > 0.5 GeV 2 
whereas for lower Q 2 values the fine DWF data appear to be higher than the data from the other 
two lattices. The solid line is a one-parameter fit form to the fine DWF data 

(Q 2 /m 2 + 1) 



GWa(Q ) = K 



(Q 2 /m 2 A + l) 2 (Q 2 /m 2 + l) , 



(50) 



which is expected assuming the validity of Eq. ( 36 ) . The fit parameter K provides an estimate of 
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the strong coupling g^NA at Q 2 = 0. In addition, we fit to the ansatz 



GWa(Q 2 



a'\ 1 - A' - 



(51) 



shown by the dashed line. The fit parameters are provided in Table III As can be seen, despite 
the fact that both fits describe sufficiently well the data for 0.5 < Q 2 < 1.5 GeV 2 , they yield quite 
different values at Q 2 = prohibiting a reliable evaluation of g^NA ■ Clearly, in order to achieve this 
goal, a better understanding of the behavior at low-Q 2 is required, since this quantity is sensitive 
to pion loop effects that maybe affected by lattice artifacts such as the finite-Ls extent. 
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blue line is a fit to pion pole dominance form of Eq. (50) for the fine DWF ensemble. The dashed black line 
is the linear fit given by Eq. (51). The strong coupling constant g n NA is the value of G^na at Q 2 = 0. 



Finally, from our lattice results we can predict the currently unmeasured ratio C^/C^, which 
is an important first approximation to the parity violating asymmetry. Its dependence in Q 2 is 



depicted in Fig. 10 From the plot we can see a very good agreement between the coarse and fine 
DWF data, at least in the range up to Q 2 ~ 1.0 GeV, indicating that there are no lattice cut-off 
effects regarding this quantity. It is also evident from the plot that at Q 2 = the ratio is expected 
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to have a non-zero value. It is noted that C]f is computed from the relationship 

r v 3 m A (m N + m A ) (r „ , , , 

and is therefore dominated by Gmi- As both C§ and Gmi lack chiral effects near the origin, the 
ratio C^jC\ is expected to be less sensitive to such effects. The present results for C^/C^ are 
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FIG. 10: The plot describes the Q 2 -dependence of the ratio C§ jC\ . The results shown are those extracted 
from both DWF lattices considered in this work. 

also consistent within statistics with the results reported earlier in Ref. |27j . 

V. CONCLUSIONS 

The nucleon to A electromagnetic, axial and pseudoscalar transition form factors are calculated 
using Nf = 2 + 1 dynamical domain wall fermions for pion masses of 330 MeV and 297 MeV for Q 2 
values up to about 2 GeV 2 . There is qualitative agreement between results obtained in the unitary 
theory and corresponding results obtained using valence domain wall quarks on a staggered sea. 
The momentum dependence of the dominant magnetic dipole, Gmi, and axial, C§, form factors 
are well described by dipole forms. They both show a slower fall-off with Q 2 than the comparison 
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to the experimental data, a fact that is reflected in the heavier dipole masses that fit the lattice 
data. Pion cloud effects are expected to dominate the low-Q 2 dependence and therefore simulations 
with pion mass below 300 MeV are required in order to allow the evaluation of such effects from 
first principles. 

The phenomenologically interesting sub-dominant electromagnetic quadrupole form factors G E2 
and Gc2 have been calculated in the case of the fine DWF lattice using the coherent sink technique 
in order to increase the statistical accuracy. The results confirm a non-zero value at low Q 2 < 
1 GeV 2 . The EMR and CMR ratios are almost Q 2 independent. The EMR values are in agreement 
with the experiment, whereas the strength of the CMR is underestimated. This can be understood 
in chiral effective theory, which predicts different chiral behavior for the two quantities. The non- 
zero values calculated in QCD are in accord with the experimental determinations [HI [381 HDTl43j 
and confirm a deviation from spherical symmetry in the Nucleon-A system. 

The axial transition form factor C§ is dominated by chiral symmetry breaking dynamics, which 
is directly reflected in the pion pole dominance. In addition, the pseudoscalar form factor G^ata 
is computed and the non-diagonal Goldberger-Treiman relation, which is a direct consequence of 
PCAC is shown to be well satisfied by the lattice data, especially for the lowest mass on the 
fine DWF ensemble. Pure monopole dependence of the Cq/C§ ratio is well satisfied, but with 
monopole masses considerably heavier than the corresponding lattice pion masses. The low-Q 2 
dependence of G w na appears to be non-trivial and the extraction of the phenomenological strong 
7r — N — A coupling, g n NA, requires careful understanding of the matrix element systematics, since 
it will be sensitive to both chiral and lattice cutoff effects. 

In conclusion, the N — A transition yields valuable information that is complementary to 
nucleon and Delta form factors. Also, since the transition is isovector, it provides an opportunity 
to assess the importance of disconnected quark loop effects. Furthermore, it provides constraints 
on the low energy constants that enter the chiral effective description of hadron properties. This 
work, utilizing dynamical chiral fermions corresponding to pion masses of 297 MeV and 330 MeV, 
together with related calculations of nucleon and Delta form factors, is a significant advance in 
the quest to understand from first principles how the closely related structure of the nucleon and 
Delta arise from QCD. The outstanding challenge for the future is to extend these calculations to 
the physical pion mass and reduce statistical and systematic errors to the level of a few percent. 
It is an appealing challenge for Lattice QCD to perform precise calculations for pion masses that 
approach the physical point with all systematics under control. Simulations with pions almost at 
its physical value will soon become available and it will be important to continue the investigation 
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of these quantities. 
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VI. APPENDIX 



Q 2 (GeV 2 ) Gmi 

DWF (N f = 2 + 1), a" 1 = 1.73 GeV, = 330 MeV 



1/11 
141 


1 ^Q1 { A(Y\ 
l.dol^4U ) 


ooU 


1. Li)o[6Z ) 


605 


0.933(33) 


819 


0.786(39) 


022 


0.641(30) 


217 


0.545(33) 


584 


0.449(50) 


757 


0.369(42) 


925 


0.332(51) 


088 


0.238(48) 


247 


0.204(99) 



TABLE IV: Coarse DWF results for Gmi, their Q 2 -dependence and the corresponding (form factor) jackknifc 
statistical errors. 



Q 2 (GeV 2 ) 


Gmi 


Ge2 


EMR (%) G C 2 


CMR (%) 




DWF (N f = 


= 2 + 1), a- 1 


= 2.34 GeV, m n = 297 MeV 


0.154 


1.602(93) 


0.0508(344) 


-3.118(2.064) 0.249(142) 


-2.748(1.595) 


0.398 


1.168(75) 


0.0146(208) 


-1.129(1.686) 0.122(98) 


-2.624(2.144) 


0.627 


0.928(84) 


0.0156(259) 


-1.528(2.749) 0.006(124) 


-3.145(4.036) 


0.844 


0.875(101) 


0.0441(375) 


-5.246(4.259) 0.158(105) 


-6.439(4.348) 


1.051 


0.593(72) 


0.0261(225) 


4.263(3.707) 0.186(67) 


-12.490(4.742) 


1.248 


0.417(86) 


0.0206(251) 


4.874(5.781) 0.197(74) 


-20.847(8.769) 


1.620 


0.439(44) 








1.802 


0.224(159) 








1.964 


0.165(181) 









TABLE V: DWF results for Gmi, G E2 , EMR (%), G C2 and CMR (%) along with their Q 2 -dependence 
shown in the first column. The errors shown are statistical jackknife errors. 
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Q 2 (GeV 2 ) 


ci 




GnNA 


DWF (Nf 


= 2 + 1), a 


1 = 1.73 GeV 


= 330 MeV 


0.141 


0.849(19) 


2.831(106) 


12.446(371) 


0.380 


0.754(19) 


1.547(55) 


13.379(436) 


0.605 


0.608(24) 


0.888(51) 


13.187(633) 


0.819 


0.604(27) 


0.755(47) 


10.941(86) 


1.022 


0.500(23) 


0.528(29) 


9.943(727) 


1.217 


0.415(26) 


0.383(28) 


9.026(870) 


1.584 


0.399(44) 


0.287(38) 


5.379(1.448) 


1.757 


0.289(38) 


0.193(28) 


5.741(1.625) 


1.925 


0.263(45) 


0.169(33) 


7.055(1.359) 


2.247 


0.186(46) 


0.093(29) 


4.942(1.740) 


DWF (Nf 


= 2 + 1), o- 


1 = 2.34 GeV 


m T = 297 MeV 


0.154 


0.825(42) 


3.103(270) 


15.292(1.005) 


0.398 


0.764(46) 


1.680(138) 


15.601(1.145) 


0.627 


0.601(61) 


0.945(139) 


11.803(1.784) 


0.844 


0.669(72) 


0.907(127) 


12.672(2.807) 


1.051 


0.502(158) 


0.579(80) 


11.556(2.098) 


1.248 


0.472(76) 


0.501(85) 


4.040(2.504) 


1.620 


0.134(278) 


-0.008(213) 


16.132(17.140) 


1.802 


0.208(161) 


0.105(117) 


3.924(6.088) 


1.964 


0.087(163) 


0.022(114) 


2.302(6.978) 


2.128 


0.097(384) 


0.084(275) 


-1.382(12.666) 



TABLE VI: DWF results for C^, C$ and G„na along with their <2 2 -dependence shown in the first column. 
The errors quoted are jackknife statistical errors. 



